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We use the fidelity approach to quantum critical points to study the zero temperature phase 
diagram of the one-dimensional Hubbard model. Using a variety of analytical and numerical tech- 
niques, we analyze the fidelity metric in various regions of the phase diagram, with particular care 
to the critical points. Specifically we show that close to the Mott transition, taking place at on-site 
repulsion [/ = and electron density n — 1, the fidelity metric satisfies an hyper-scaling form which 
we calculate. This implies that in general, as one approaches the critical point U = 0, n — 1, the 
fidelity metric tends to a limit which depends on the path of approach. At half filling, the fidelity 
metric is expected to diverge as [7~* when U is sent to zero. 
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I. INTRODUCTION 

Recently a novel characterization of phase transitions 
has been advocatecP^. This is the so called fide lity ap- 
proach to critical phenomenePE^^ 5 6 7 8 9 lo ii 12 I3 ii4|i5 | i6 | 

that relies solely on the state of the system and does not 
require the knowledge of the model Hamiltonian and its 
symmetry breaking mechanism. Two states of the sys- 
tem at nearby points in parameter space are compared 
by computing their overlap (the fidelity). Since quantum 
phase transitions are major changes in the structure of 
the ground state, it is natural to expect that, when one 
crosses a transition point the fidelity will drop abruptly. 
To make the analysis more quantitative one considers 
the second derivative of the fidelity with respect to the 
displacement in parameter space. Remarkably this sec- 
ond derivative, more in general the Hessian matrix, de- 
fines a metric tensor (the fidelity metric hereafter) in the 
space of pure state^^^ A super-extensive scaling of the fi- 
delity metric corresponds to the intuitive idea of a fidelity 
drop. Indeed, it was shown that, at regular points 
the fidelity metric scales extensively with the system size, 
and a super-extensive behavior implies criticality. How- 
ever the converse is not true in general; in order to ob- 
serve a divergence in the fidelity metric a sufficiently rel- 
evant perturbation (in the renormalization group sense) 
is needecP^. Loosely speaking the more relevant the op- 
erator the stronger the divergence of the fidelity metric. 
The Berezinskii-Kosterlitz-Thouless (BKT) transition is 
peculiar in this respect as it is driven by a marginally rel- 
evant perturbation, i.e. with the smallest possible scaling 
dimension capable of driving a transition. This gives rise 
to an infinite order transition and as such the BKT does 
not rigorously fit in the simple scaling argument given 

iijni. 

Surprisingly, contrary to the nai've expectation, Yang 
has showijISthat in the particular instance of BKT tran- 
sition provided by the XXZ model, the fidelity metric di- 



verges algebraically as a function of the anisotropy. This 
is an appealing feature since observing a singularity at a 
BKT transition is generally a difficult task- 

In this paper we analyzed, with a variety of analytical 
and numerical techniques, the one dimensional Hubbard 
model primarily aiming at assessing the power and lim- 
itations of the fidelity approach for infinite order QPT 
(n = 1, [/ ^ 0). We believe it is useful to list here the 
main accomplishment of our analysis: i) An exact cal- 
culation, on the free gas line U — 0, shows that the 
fidelity metric g presents a cusp at half filling and a 
1/n divergence at low density n respectively, ii) Using 
Bosonization techniques, we observe a divergence of the 
form g in the regime U n when n ^ 0. iii) 

Resorting to Bethe Ansatz we are able to interpolate 
between the regime where the Luttinger liquid param- 
eter K approaches the BKT value 1/2 and that where 
K ^ 1, which describes the free-Dirac point. We show 
that the fidelity metric satisfies an hyper-scaling equation 
which can also be extended to finite sizes, iv) We cal- 
culate the hyper-scaling function in the thermodynamic 
limit by solving Bethe-Ansatz integral equations while 
at half filling by resorting to exact diagonalization. As 
a consequence, when approaching the transition point 
U = 0, n — 1 the fidelity metric tends to a limit which 
depends on the path of approach. On the particular 
path U ^ 0, n — 1, an algebraic divergence of the form 
g ^ U~'^ is expected, on the basis of numerical results. 

In the ID Hubbard model the BKT transitiorP^ occurs 
exactly at half filling as soon as the on site interaction U 
is switched on, inducing a gap in the charge excitation 
spectrum. Away from half filling instead all modes are 
gapless for any U and the system is a Luttinger liquid. 
Since at half filling, the only gapless point is at J7 = 0, 
the kind of BKT transition offered by the Hubbard model 
is different from that featured by the XXZ model. In 
that case one continuously arrives at the transition point 
from a gapless phase by tuning the anisotropy parameter. 
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This difference, in turn, makes more difficult the analysis 
of the fidelity metric in the Hubbard model. 



II. PRELIMINARIES 



The one dimensional Hubbard model is given by 



H = -t 



E (4. 



Ci+i,CT +h.c.j +[/ ^ 



(1) 

where n^.o- = cl ^Ci^cr , n-i = ni.| + n^.j and we will be 
concerned with the repulsive / free gas case when U >0. 
Due to the symmetries of the modeP^ it is sufficient to 
limit the analysis to filling small er or equal to one half 
i.e. n = (rii) < 1. It is well knowrPsHSl that for n < 1 the 
model is in the Luttinger liquid universality class for any 
value of the interaction U. Spin and charge degrees of 
freedom separate and their respective excitations travel 
at distinct velocities Vg and Vc- Both charge and spin 
modes are therefore gapless. Exactly at half filling (n = 
1) the system becomes an insulator and develops a charge 
gap AE^ ^fi+-^i_^E{N+l)~2E{N) + E{N- 1), 
where E (N) is the ground state energy with N particles. 
The gap opens up exponentially slow from U = 0, and the 
point U — 0, n — 1 is a transition of BKT typ^2l]_ The 
length scale ^ = 2vc/ AE^, describes the size of soliton- 
antisoliton pairs, in the insulator. As we approach the 
critical point at C7 = 0, n = 1, these pairs unbind and 
proliferate, allowing the system to conduct. 

In the fidelity approach one is interested in the 
overlap between ground states at neighboring points 
of the coupling constants (say a vector A): F {X) — 
(A) {ip (A + dX))\. Remarkably, the second order term 
in the expansion of the fidelity defines a metric in the 
space of (pure) states: 



F(A) = 1 



1 



where 
G 



Re [{df,tpQ\di,ipo) - {df,ipo\tpQ){tpQ\di,ipo)] 



and tpQ = ip (A) and = d/dX'^. Actually, at regu- 
lar points A of the phase diagram, G^,!^ is an extensive 
quantitj^H so that it is useful to define the related in- 
tensive metric tensor ,y = G^,u/L. With reference to 
Hamiltonian ([ij it is natural to investigate the behavior 
of the fidelity under variations of the interaction param- 
eter U . The possibility of analyzing variations of the 
chemical potential, though appealing does not fit in the 
fidelity approach as the ground states at different be- 
long to different super-selection sectors. Hence now on 
we will solely be interested in guu and we will simply 
write g in place of gjj.u- In Ref .II^] it was shown that it 
can be written in the following form 



Wn 



V 



(2) 



where En, \n) are the eigen-energies and corresponding 
eigenstates of the Hamiltonian ([iJ (with repulsion U and 
filling n), |0) corresponds then to the ground state and 
Vi,j — {i\V\j). In passing, we would like to notice that de- 
spite the apparent similarity between Eq. ^ and the sec- 
ond derivative of the energy E" (A) (a similarity stressed 
in Ref.ISS])^ it is possible to show - using the Rayleigh- 
Schrodinger series - that the metric tensor is in fact re- 
lated to the third (and first) energy derivative, via 



G = 



1 



Vbo ^ (Er 

i,j>0 ^ 



E 



Vo.iVijVj.O 
— Eq) {Ej — Eq) 



E'" 

1^' 



Moreover, in cases where E'" (A) is bounded 
in the thermodynamic limit, one obtains 
the interesting kind of factorization relation 

{o\v\o){o\vg{Eofv\o) = {o\vg{Eo)vg{Eo)v\o) 

valid in the thermodynamic limit, where Q is the 
resolvent g {E)={1 - \0){0\) {H - E)~'^ {1 - |0)(0|). 

In the rest of the paper we will be concerned with the 
analysis of the metric tensor with special care to the BKT 
transition point. The phase diagram of the Hubbard 
model is depicted in Figjl] The model has been solved 
by Bethe Ansatz in Ref.lSST We will tackle the problem 
using a variety of techniques. On the free-gas U = line, 
an explicit calculation is possible at all fillings. Around 
the region U = and filling away from n = and n = 1 
bosonization results apply. Instead, close to the points 
U — and n = 1 we will cross results from bosoniza- 
tion with Bethe- Ansatz in order to extend bosonization 
results up to the transition points. We will show that the 
behavior of the metric is encoded in a scaling function. 
Away from half filling the scaling function is computed 
integrating Bethe-Ansatz equation, while at half filling 
by resorting to exact diagonalization. 



III. EXACT ANALYSIS AT U = 

At U — the complete set of eigenfunctions of Hamil- 
tonian ([iJ is given by a filled Fermi sea and particle-hole 
excitations above it. It is then possible to apply directly 
Eq. (§. 



A. Half filling 

We first treat the half filled case n = 1, 
where the Fermi momentum lies at kp = tt/2. 
Writing the interaction in Fourier space as y = 
c^_|_^ tCfc,tCfe',j.,and going to the ther- 



modynamic limit, Eq. p.\ becomes 



_ 1 f rik (1 - rtfc+g) rik' (1 - Uk'- 

(27r)^ J[~Tr,Trf (ek+q - Efc + Efc'-q - Cfe' ) 



(3) 
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Figure 1: Phase diagram of the repulsive Hubbard model. 
The region < n < 1 is in the Luttinger liquid (LL) uni- 
versality class. The hatched area corresponds to the condi- 
tion 5^ iU) < 1, where ^ (U) is the correlation length at half 
filling. Approaching the critical point (n = 1,(7 = 0) from 
within (the complementary of) this region the LL parameter 
Kc approaches 1/2, the BKT value (1, the free-Dirac value). 



where Ck — — 2tcos(fc) is the U — single par- 
ticle dispersion and = e^) {— Uk^f — f^A:,x 
in absence of magnetic field) are the fermionic, zero- 
temperature, filling factors. Using e{k,q) = €k+q — ek — 
4sin {q/2) sin (fc + q/2) and substituting k' —k' we ob- 
tain 



dq- 



1 



(27r)^7o 8sin(g/2)^ ..-^ 



dk 



dk' 



n{k,q)n{k',q) 



[sm{k + q/2)+sm{k' + q/2)Y 



Changing variables to p = k + q/2, p' = k' + q/2 and 
making a shift of 7r/2 we obtain finally 



(27r)^yo 8sin(g/2)^ 



1 



247r2 



0.00422 (4) 



where we defined 



Jiq) = 4 



q/2 



q/2 

dp / dp' - 



1 



JO [COS (p) -t- cos (p')]" 
In (cos (q/2)) 
cos (q) — 1 ' 



and correctly J (q) = J {—q) > 0. 

A related interesting issue is that of the finite size 
scaling of the metric tensor g i.e. the way in which 
at length L converges to its thermodynamic value (see 
alscPl). In Ref.li^ it was shown that in a gapless regime 
scaling analysis predicts g^ ~ apart from regular 

contributions which scale extensively (and contribute to 
g with a constant). Here = 2Av — 2C — 1 where Ay is 
the scaling dimension of V in the renormalization group 



sense and is the dynamical critical exponent. On the 
line U = one has Ay = 2 as F is a product of two 
independent free fields, while ^ — 1 when n ^ due 
to the linear dispersion of excitations at low momenta. 
This implies that at leading order g^ ^ A + BL~^ . One 
should however be careful that logarithmic corrections 
are not captured by the scaling analysis of Ref.^i^ and 
they might be present due to the BKT transition occur- 
ring at this point. Let us try to clarify this issue. Looking 
at Eq. ([4]), as a first approximation, we might think that 
the finite size g^ is well approximated by the Riemann 
sum of F{q) = J (q) / sin {q/2f: 



9L 



Sl 



(27r) 



St = 



2tt 



L/2 



L ^ 

?i=i 



F 



2tt 



Now F (q) diverges logarithmically around tt : F (q) 



-4 In I 



O 



and since the Riemann 



sum of the logarithm converges to its integral as 
{A + BlnL) /L + (L-^) we would conclude 



9L 



1 

247r2 



A- 



BlnL 



C 
1? 



0{L- 



(5) 



However, a detailed analysis of the finite size g^ reveals 
that a cancellation occurs between two logarithmic cor- 
rections so that actually S = in Eq. ([5j. The exact 
finite size g^ is composed of two terms. One term is 
a triple sum which, in the thermodynamic limit, corre- 
sponds to the triple integral in Eq. ([3|. The other term is 
a double sum which originates from zero transferred mo- 
mentum contribution and vanishes when L — > 00. We nu- 
merically verified that both terms contain a In L/L part 
when L 00, but their contribution is equal and oppo- 
site so as to cancel out exactly from g^. The absence of 
logarithmic corrections can clearly be seen in the inset of 
Fig. [2] where the finite size is plotted against 1/L. 



B. Away from half filling 

Similar considerations can be done away from half- 
filling. Eq. ([3} still holds, simply in this case kp ^ 7r/2. 
We assume kp — n7r/2 < tt/2 (anyway a particle hole 
transformation implies g {n) = g{2 — n)). First note 
that the integral over q can be recast as 2 JJ^ dq. The 
filling factors constrain the momenta to |fc| < kp and 
|fc -I- g| > kp. If q < 2kp this implies kp — q < k < kp. 
Instead if q > 2kp the sum over k is unconstrained: 
—kp < k < kp. Thus we obtain 

2 1 n'2kp' pkf rkp 

^ = 7^ 1 / '^'^ / 

(ZTT) yjo Jkp-q Jkp-q 

/kp pkp I 
dk I dk' \ 
^ -kp J -—kp J 

1 

16 sin (q/2)^ [sin (fc + q/2) + sin (fc' + q/2)]^ ' 
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Changing variables as before and defining 



J (a, 6) 
we obtain 

9 = 



nb pb 

I dp / dp'- 

a J a 



1 



+sin {p')Y 




J{kF-q/2,kF + q/2) 



sm{q/2y 
J{-kF + q/2,kF 



9/2) 



sin(g/2)^ J 

In Fig. [2] one can see a plot of g{n,U — 0) as a function 
of the total density n — Ntot/L = 2fcp'/7r. It is possible 
to show that in the very dilute regime n ^ 0, the fidelity 
metric g diverges in a simple algebraic way 

5(n^0,C/ = 0)~-. 

n 

This divergence can also be simply understood by re- 
sorting to the scaling arguments reported irP3] There 
it was shown that, in the thermodynamic limit g ~ 
~ Mcl^"^^^ where now is the scaling dimension 
of the field jj.. On the line J7 = as already noticed 
Ay = 2 while now C = 2 as n ^ to account for the 
parabolic dispersion. The chemical potential scaling ex- 
ponent is A^ = 2 in the dilute Fermi ga^^. All in all 

Ui • I 1-1/2 -1 • I |l/2 

we obtam g \^ — ^ n smce n ~ |/i — /id 

which agrees with the explicit calculation. 

The finite size scaling of g^ for different filling < n < 
1 is the same as that observed at n = 1 and is dictated 
by Eq. ([s]) with i? = 0, as can bee seen in the inset of 
Fig.|2] 

Finally note that, since g {n,U = 0) is symmetric 
around n = 1, g (n) has a local maximum at that point 
with a cusp. The origin of this discontinuity is not well 
understood at the moment but reveals a signature of the 
transition occurring at this point. 



IV. BOSONIZATION APPROACH 



It is well knowrJ^I ^^ I ^^ I that for U > and away from 
half filling (n = 1), the low energy, large distance behav- 
ior of the Hubbard model, up to irrelevant operators, is 
described by the Hamiltonian 



H = H, 



d^3 



iv,. 



(6) 
s, c . 



Charge and spin degrees of freedom factorize and are de- 
scribed respectively by He and Hg. The Luttinger liquid 
parameters K^^s are related to the long distance, alge- 
braic decay of correlation functions, while Vc.s are the 
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Figure 2: Fidelity metric g as a function of the total density 
at U = 0. The singularity at n — > is of the form . In the 
inset the finite size scaling oi Ql for some different fillings is 
shown. The approach to the thermodynamic value is given by 
Eq. ([s]) with 5 = 0. In fact fitting the data points with Eq. ^ 
gives values of B/A of the order of 10~^ and a chi-square of 
the same order as the one obtained with B = 0. 



speed of elementary (gapless) charge and spin excita- 
tions. From bosonization, and setting the lattice constant 
a = 1, one finds, for small U, 



Kr 



1 



U 



2'KVl 



(7) 



where the Fermi velocity is vp = 2tsin(fci?) and kp = 
■Kn/2. Instead the Luttinger parameter Kg is fixed to 
Ks = 1 due to spin rotation invariance. Exactly at n = 1 
there appears another term (an Umklapp term) in the 
charge sector which is marginally relevant and is respon- 
sible for the opening of a mass gap. In this case the 
effective theory is the sine-Gordon model. 

Since the fidelity of two independent theories factorizes 
the metric tensor g is additive and we obtain g = gs +gc- 
In Ref.^ the fidelity metric of a free boson theory has 
been calculated to be given by 



1 / 1 dK,, 



8 dU 



(8) 



In our case = as Kg does not vary so that g = 
gc + gs = gc- Using Eq. ([tJ one obtains a formula valid 
up to zeroth order in U : 



1 



2{4:TtvfY 



0{U). 



(9) 



Some comments are in order. The expansion ([tJ is ac- 
tually an expansion around J7 = valid when U ^ vf- 
When we move toward the BKT critical point one has 
^ 2 and g ^ 1/(12877^). This value is 3/16 the 
number calculated directly at J7 = in Sec. [Ill A[ We 
believe that this discrepancy is due to lattice corrections 
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which are neglected in formula ([s]). Approaching the low 
density critical point U = 0, n = Eq. ([9} predicts that, 
in a narrow region U ^ n, the fidelity metric diverges 
as (7 ~ n^'^. This contrasts with the result g ^ 
obtained at U — 0, as one would expect since U is a 
relevant perturbation. In fact, in the diluted regime, 
the low-energy effective theory is that of a spinful non- 
relativistic gas with delta interaction^^. There one still 
has n ~ l/i — fJ.c\^^^ and dynamical exponent C = 2. Thus 
if we take — 1, then using the conventional scal- 
ing analysis we would find g ~ n~^, consistent with the 
bosonization result. 



V. HYPER-SCALING OF FIDELITY METRIC 
NEAR THE METAL-INSULATOR CRITICAL 
POINT 

The bosonization expression Eq. ([?]) is an expansion of 
Kc {n, U) around U = where Kc reaches its free Dirac 
value 1. In the whole stripe U > 0, < n < 1, Kc (n, U) 
is a bounded function ranging between 1/2 and 1123 The 
maximal value Kc = 1 is obtained in the segment U — 0. 
Instead the minimal value Kc = 1/2 is attained at the 
lines n = and n = 1 and in the strong coupling limit, 
i.e. Kc ^ 1/2 for U ^ \t\. This considerations show that, 
from Eq. (|8|, g can be infinite only at the points U = 
and n ~ 0, or 1 where Kc is discontinuous. In particular 
we are interested to the vicinity of the transition point 
U = 0, n — 1 which we will call simply (with some abuse) 
BKT point. Calling ^ (U) the correlation length at half 
filling, and 6 — 1 — n the doping concentration, it can be 
shown (see later) that the Luttinger liquid parameter Kc 
tends to 1/2 when the BKT point is approached from the 
region 5^ (U) <C 1. Instead Kc 1 when the BKT point 
is approached from S^{U) ^ 1. Given this discontinuity 
of Kc it seems difficult to interpolate between the two 
regimes S£, (U) ^ 1 and 5£, (U) ^ 1. However we will 
show that such an interpolation is indeed possible and 
that the fidelity metric g satisfies an hyper-scaling rela- 
tion valid when (5^ ranges over many order of magnitudes 
(see figure [3|. 

In Refs .122122] g^jj efficient characterization of the Lut- 
tinger liquid parameter Kc in terms of Bethe Ansatz 
results has been found. Kc is related to the so called 
dressed charge function Zk through 



Kc 



4/2. 



(10) 



where the wave vector Q is a generalization of the Fermi 
wave- vector in the interacting regime, and has to be de- 
termined by Bethe- Ansatz equations. We will now argue 
that, in the metallic phase, the dressed charge function 
Zq satisfies the following scaling relation: 



ZQ{U,S)^<i>z{HU) S). 



(11) 



Here ^ (U) is the correlation length at half filling {6 = 0) 
defined via ^ = Vc/AEc where Vc is the charge carries 



(holons) velocity and AEc is the (charge) gap, and 
is a scaling function. Obviously a similar scaling relation 
holds for the Luttinger parameter Kc through Eq. (10 1. 



The scaling relation holds as long as the interaction is 
not too strong, say U < 1. A similar scaling relation has 
been conjectured in Ref.'22] for the charge stiffness. 

Using equation (|8} we obtain the following scaling re- 
lation for the fidelity metric of the Hubbard model in the 
metallic phase close to the BKT point 



9iU,6) = 



dine 
dU 



(12) 



where we introduced the scaling function <i>g (x) ~ 
(1/2) {x<^'z i^) I'^z {x)f. Following Ref.l2al it is natural 
to conjecture a more general hyper-scaling relation for 
the fidelity metric valid also at finite size. Building the 
other dimensionless quantity with the correlation length 
e and the size L, we obtain 



g{U,5,L) 



dU 



(13) 
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Figure 3: (color online). Scaling function <I>g for the fidelity 
metric as a function of the scaling parameter ^5. Data are 
obtained by solving numerically the integral equations for the 
dressed charge at different interaction strength U. The dashed 
line is obtained expanding the Bethe- Ansatz equations for the 
Luttinger parameter in the regime ^5 ^ 1 up to second order 
in 5. The solid line results from integrating the RG BKT 
equations and is valid when ^5 ^ 1. 



We will now verify the hyper-scaling relation Eq. ([13| 
away from half-filling in the thermodynamic limit {£^/L — 
0) using various analytical techniques and at finite size 
resorting to exact, Lanczos, diagonalization. 



A. Away from half filling 



The scaling relation Eq. ( [T2| (and implicitly Eq.(lll) 
can be verified analytically in the two limits x = £,6 ^ 
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and X oo. We present the results in terms of the Lut- 
tinger param eter Kc which has more physical relevance. 
In RefsP^^, the Bethe Ansatz equations for the dressed 
charge have been solved around ^ = 0. At leading order 
in 6, they found for the Luttinger parameter 



K,^l + aiU)S+l[a{U) 6f + O {S^) 



(14) 



The function a ([/) is studied in the appendix and is 
approximately given by the following expansion, for 
U/2t: < 1 



aiU) 



In (2) 



,27v/U 



u 



Not surprisingly, this has the same form as the soliton 
length, ^ ([/). In fact, for the regime of interest, U/27r < 
1, Vc ^ 2 and ^{U) = na{U) /21n(2j^. This implies 
that in the region x ^ 1, the metric scaling function 
behaves, at leading order, as 



1>9 i^) = 



In 4 
V2Tr^ 



0\ 



Conversely, in the opposite regime ^5^1 (i.e. U 
0, S small) integrating the renormalization-group BKT 
equationP^, one is able to improve the bosonization re- 
sult Eq. Q and we obtain 



U/ (47r) 



1 



- C//(27r) ln(l/7r5) 
[ln(7re<5)]-V.... 



Accordingly, the scaling function <i>g (x) has the following 
asymptotic behavior when a; ^ oo 



(^) 



1 

32 



[In (ttx)] 



In the limit J7 ^ we recover bosonization's result g 
1/I287r2. 

To verify the scaling relation Eq. (12 I also in the inter- 
mediate regime ^5 « 1, we solved numerically the Bethe- 
Ansatz equations for the dressed charg^^^. To obtain the 
dressed charge function we need also the density of 
wave numbers pk- They are solutions of the following 
integral equations 



f 

Pfc = 1 — cos (fc) / dq COS qR {sin k — sin q) pq (15) 
J-Q 

rQ 

Zk = I + / dq cos qR (sin k — sin q) Zq (16) 
J-o 



with R (x) given by 
Rix) = 



1 

27r 



duj 



The wave- vector Q is determined by fixing the electronic 
density n = J^^ pkdk. Integrating numerically Eqs. (15 I 



and ( 16 1 for different value of the coupling strength and 
doping fraction, we are able to verify the scaling relation 
Eq. (121 over many order of magnitudes. The result is 
plotted in figure |3] 

We would like to point out that, since close to the BKT 
point the relevant variable is ^ (U) S, the limit of the fi- 
delity metric when U ^ 0, n ^ 1 depends on the path of 
approach. However, as we have shown, the combination 
g {d\n£_/dU) ^is a perfectly well defined function of £_5. 



B. Numerical Analysis at half filling 

To study the scaling behavior of the metric at half- 
filling, we turned to exact, Lanczos, diagonalization. We 
find that the metric obeys the scaling form Eq. ( 13 1 with 
the scaling function Y (0, y) plotted in Fig. (|4]), for values 
of its argument ranging over six order of magnitudes. 

After calculating the ground state (U) of Hamilto- 
nian ([ij with the Lanczos algorithm, the intensive fidelity 
metric g is obtained from the fidelity F{U, U + SU) — 
li^omi'^oiU + SU))\ using 



9L 



2 1 



F{U, U + SU) 
61P 



(17) 



with 6U = 10^'^. The above equation is a good ap- 
proximation to the limit SU ^ as long as SU <C 
1/^/Lg which was confirmed to be the case in all sim- 
ulations. The function Y {0,y) is then obtained via Y — 
gL{dln^/dU)~ . In order to have as many data points as 
possible, we analyzed sizes of length L = An+2 with peri- 
odic boundary conditions (BC's), while for length L = 4n 



Y(0,y) 









' 10 




' 12 




' 14 




- (l/6)ln(Cyf C^9.92 


: 





Figure 4: (color online). Finite-size scaling function Y {0,y) 
as a function of the normalized correlation length ^/L. The 
solid line is obtained by requiring consistency with the non- 
interacting value at ?7 = 0. The constant C is obtained by a 
best fit with the numerical data (symbols). Boundary con- 
ditions are antiperiodic when system size L is a multiple of 
four, periodic otherwise. 
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antiperiodic boundary conditions were used. Choosing 
such boundary conditions at half fiUing, assures that the 
ground state is non-degenerate even at [/ = Cp3 

The behavior at large y = £_/L can be obtained by 
requiring consistency with the value obtained in the free 
case U — Q. Then the scaling function must have the 
following limiting form 



y -> oo, 



where C is a constant. Since when y ^ oo In^ ~ 2-k/U 
the divergence in (dln^/dC/) cancels with the logarithm 
above, and one obtains 



9{U-^Q,L) 



1 



247r2 \ \n{C^/L) ) 
1 



247r2' 



Having computed the scaling function Y {x,y) we could 
ask what happens to the metric tensor as one approaches 
the BKT point from the particular path U ^ 0, n = 1. 
The smallest values of y at our disposal are of the order 
of y '--^ 10~^. Looking at Figure [i] on the basis of these 
data, it seems that the function Y approaches a non-zero 
value limj,^o Y {Q,y) as y goes to zero. If this is the case, 
after taking L oo and U small we would have 



lim g ([/, n — 1,L) 



r(o,o) 



Note however, that observing this divergence numerically 
can be very hard as we must be in the region L ^ ^ ([/) 
which requires huge sizes when the coupling U is small. 



VI. CONCLUSIONS 

In this article we analyzed the fidelity metric in 
the zero-temperature phase diagram of the ID Hub- 
bard model, with particular care at the phase transition 
points. The fidelity metric quantifies the degree of dis- 
tinguishablity between a state and its neighbors in the 
space of states, and as such it is expected to increase 
(or diverge) at transition points. Special attention has 
been drawn to assess whether the fidelity metric reveals 
signatures of the Mott-insulator transition occurring at 
on site repulsion f/ ^ and filling factor n — \. Being 
a transition of infinite order, it is particularly difficult 
to pin down since typical thermodynamic quantities are 
smooth (although not analytic) at the transition. The 
point [/ = 0, n = 1 is particularly singular in that it 
represents the limit of two completely different physical 
regimes. On the line [/ = it is simply the half-filling 
limit of a gapless free system, whereas fixing rt = 1 it 
represents the limit of a complicated interacting massive 
system. Surprisingly, we have shown that it is possi- 
ble to interpolate between these two regimes, and the 



fidelity metric defines a hyper-scaling function which de- 
pends only on a; = ^ {U){1 — n). The two regimes roughly 
correspond to x ^ land x ^ 1 respectively. 

Away from half filling we have been able to com- 
pute the scaling function integrating numerically Bethe- 
Ansatz equations, and we obtained analytic expressions 
for the limits x 0,oo. The result implies that, as a 
function of U and n separately, the fidelity metric has 
no precise limit when J7 — > 0, n ^ 1, but the scaling 
function is well defined in term of the scaling variable x. 

Precisely at half filling we computed the scaling func- 
tion resorting to exact diagonalization and upon intro- 
duction of another scaling variable y — {U) /L. With 
the numerical data at our disposal, the scaling function 
appears to be smooth and non-zero around y = 0. As a a 
consequence, approaching the Mott point from the half- 
filling line, the fidelity metric would display a singularity 
of the form U~*. A singularity of algebraic type has been 
observed also in another instance of Kosterlitz-Thouless 
transition given by the XXZ model. 

LCV would like to thank Cristian Degli Esposti Boschi 
for a critical reading of the manuscript. The work of 
P.B. has been partially carried out in the framework of 
the PRIN project Microscopic description of fermionic 
quantum devices. 



Appendix A: 
Following Ref.f^Sl the coefficient a in Eq. ( 14 1 is given 



by 



where 



a{U) 



4 In (2) 

uJW) 



and Jo is a Bessel function. Using the following results 
1 



)OLX 



n=0 



1 



Jo (a;) e-'^^dx = , 

v^+w 



, a,(3>0, 



we arrive at 



OO 

/(f/) = l-|-2^(-l)"^= 



Using the formalism of the Remnant Functions defined in 
Ref.l^one realizes that / ([/) is related to Ry^ ^ (4[/~2). 

With the help of the expansions in Ref.l2^ we obtain the 
following expression 



n=0 ^ 
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It is now easy to obtain the desired expression, using the 
asymptotic of the Bessel function 



Ko {II x) = 




Collecting the results together, we obtain at leading order 
a([/) = M£)e2VC/ (1 + ([/)). 
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ness which experiences a jump at the transition (see i.e.l^ 
The term BKT is justified in the sense that the underlying 
effective theory is the sine-Gordon model (see later). In 
the condensed matter literature the term Mott transition 
is more commonly used. 
^'^ Note also that both these classes of fermionic systems with 
different boundary conditions, can be mapped onto an in- 
teracting systems of hard-core bosons with periodic bound- 
ary conditions. This clarifies further why it is possible to 
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